Non-Markovian probes in ultracold gases 
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We present a detailed investigation of the dynamics of two physically different qubit models, 
dephasing under the effect of an ultracold atomic gas in a Bose-Einstein condensed (BEC) state. We 
study the robustness of each qubit probe against environmental noise; even though the two models 
appear very similar at a first glance, we demonstrate that they decohere in a strikingly different 
way. This result holds significance for studies of reservoir engineering as well as for using the qubits 
as quantum probes of the ultracold gas. For each model we study whether and when, upon suitable 
manipulation of the BEC, the dynamics of the qubit can be described by a (non-)Markovian process 
and consider the the effect of thermal fluctuations on the qubit dynamics. Finally, we provide an 
intuitive explanation for the phenomena we observe in terms of the spectral density function of the 
environment. 



I. INTRODUCTION 



All realistic quantum systems interact with an envi- 
ronment and a proper description of their dynamics calls 
for the toolbox of open quantum systems theory [l|, Q ■ 
In this framework the simplest type of noise is described 
by a Markovian master equation in the Lindblad form, 
corresponding to a completely positive, trace preserving 
dynamical map satisfying the semi-group property Q. 
The latter condition means that the map can be divided 
into infinitely many time-steps, each identical and inde- 
pendent of the past and future steps Q, and therefore a 
Markovian dynamical map has the intuitive interpreta- 
tion of memoryless dynamics. Markovian processes suc- 
cessfully describe a plethora of physical processes, partic- 
ularly in the field on quantum optics, but can fail if ap- 
plied to more complex system-environment interactions 
where memory effects become important. In such situa- 
tions one must resort to non-Markovian dynamical maps. 
Recently the theory of non-Markovian dynamics has 
beautifully taken shape as a result of proposals for the 
definition of non-Markovian dynamics [JdMl- Amaz- 
ingly, in the past the whole concept has lacked a sim- 
ple, model-independent definition. The application of 
non-Markovianity quantifiers, constructed on the basis of 
these definitions, has led to a deeper understanding of the 
microscopic mechanisms underlying non-Markovian dy- 
namics [l(| , and to the clarification of hazy concepts [ll| . 
Moreover, it turns out that the quantifiers can be used 
to witness initial correlations in the composite system- 
environment state [l2|, and in the environment state, 
and to probe quantum phase transitions of the environ- 
ment [13l | . to name just a few examples. The quantifier 
put forward by Breuer et al, equating non-Markovianity 
with bidirectional information flow, has also been stud- 
ied in a linear optics set-up, thus establishing that non- 
Markovianity quantifiers are not merely a theoretical tool 

In the spirit of these advances we recently conducted an 
investigation of the non-Markovian dynamics of a qubit 



coupled to an ultracold Bose-Einstein condensed (BEC) 
gas with a two-fold motivation [i"o| . On one hand, exper- 
imentalists have discovered astonishingly accurate means 
of controlling and manipulating ultracold gases ■ This 
raises a question whether ultracold gases could provide a 
tailored environment to a quantum system such that its 
decohering effect on the system is minimized. Indeed, we 
discovered that simple and experimentally feasible ma- 
nipulation of the ultracold reservoir leads to significant 
changes in the way the qubit dephases and enables a 
perfect control of the Markovian to non-Markovian tran- 
sition in the qubit dynamics. On the other hand, the 
way a qubit decoheres may reveal important information 
about the environment, leading to the concept of a probe 
qubit [l7l - [l9| . Here the fundamental question is to what 
extent one may probe a large, complicated environment 
by looking at a simple and accessible quantum system 
that interacts with it. This work aims to dive deeper 
into these two aspects in the context of qubits embedded 
in ultracold gases. 

More specifically, we study the dynamics of two differ- 
ent qubit models, each embedded in an identical envi- 
ronment, namely a BEC gas. The scattering length of 
the bosons forming the BEC can be controlled using the 
Feshbach resonances and therefore we have access to an 
environment that can be chosen to consist of cither free 
or interacting bosonic particles. It is worth stressing that 
the latter regime is widely unexplored and prototypes of 
open system models are mainly built on the assumption 
of an environment of free particles. Intuitively one may 
expect the interacting environment to have better mem- 
ory keeping properties than a non-interactin g on e and 
this was exactly what we discovered in Ref. [15| : non- 
Markovian effects take place when the environment is 
sufficiently strongly interacting. It was left as an open 
question, however, whether this phenomenon is specific 
to the model we studied in Ref. |l5| or if one can gener- 
ally associate interacting environments to non-Markovian 
dynamics. By comparing and contrasting the reduced dy- 
namics of two different qubit models we find the answer 
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to be negative: an interacting environment can induce 
Markovian dynamics on one qubit architecture and non- 
Markovian on another. 

We also address another unresolved question of non- 
Markovian open quantum systems, namely the connec- 
tion between the emergence of memory effects and the 
form of the spectral density function characterizing the 
dynamical map. Non-Markovianity is often associated to 
structured spectra. In the case of the Jaynes-Cummings 
model, for example, a decrease in the width of the 
Lorentzian spectral density function always leads to a 
higher degree of non-Markovianity [20|. The connection 
is much more subtle for purely dephasing processes, as 
shown in Ref. pH | , where two of us unveil a condition on 
the form of the spectrum to create non-Markovian dy- 
namics. In this article we study this connection for the 
rather complex dephasing dynamics induced by the BEC 
environment and show that, quite unexpectedly, purely 
Markovian dynamics can arise from spectral density func- 
tions with rich structure. 



II. THE TWO MODELS 
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Model I: Double well qubit Model II: Atomic quantum dot 

FIG. 1. (Color online) The two qubit architectures considered 
in this article. Model I comprises of an impurity atom trapped 
in a double well potential and model II assumes a single impu- 
rity atom, with an internal level structure, trapped in a deep 
harmonic trap. Each qubit interacts with a Bose-Einstein 
condensed ultracold atomic gas in a shallow harmonic trap. 



We compare two different qubit models composed of 
a trapped impurity atom interacting with an ultracold 
bosonic gas. In model I, originally introduced in Ref. [22l | 
and displayed in Fig. QJa), the impurity atom is trapped 
in a deep double well potential and forms an effective 
qubit system where the two qubit states are represented 
by occupation of the impurity in the left \l) or the right 
well |r). In model II, see Refs. pHH and Fig. [TJb), the 
impurity is trapped in one site of an optical lattice and 
has two internal states, |e) and \g), representing the qubit 
states. The Hamiltonians for both models are composed 
of three parts: Hamiltonians of the impurity and of the 
interacting background gas and the interaction Hamilto- 
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nian, respectively, 

/2 
dx*t( x ) J^L + Vb(x) + ^$t (x)$(x) 

Hab = ~7r~ [ dx$t(x)vKt( x )v]/(x)$(x). 



<f(x), 

(1) 



Here tua, ^( x ) an d Va(x) ar< 3 the mass, field opera- 
tor and the trapping potential of the impurity atom, 
nig, $(x), gg = A"K% as/mB and Vb( x ) ar e the mass, 
field operator, coupling constant and the trapping po- 
tential of a background gas atom and as is the scatter- 
ing length of the boson-boson collisions. Finally, g A B = 
4Trh 2 aAB /mAB is the coupling constant of the impurity- 
boson interaction where m A B — rnArriB / '(jua + tub) is 
the effective mass. 

In both models we expand the impurity held operator in 
terms of Wannier functions {</>k} localized in the lattice 
sites/the two wells. Assuming that the lattice sites/the 
two wells are very deep, hopping and tunneling effects 
are both suppressed and the Wannier functions take a 
Gaussian form. We assume that the background gas is 
weakly interacting and can be treated in the Bogoliubov 
approximation, neglect all terms that are quadratic in 
the creation and annihilation operators of the Bogoliubov 
modes and assume that the background has is homoge- 
nous. 

It turns out that when we focus on a single impurity 
the Hamiltonians Ha and Hg in moth models arc effec- 
tively the same. Any differences in these Hamiltonians 
will not have an effect on the dynamics of information 
flow characterising non-Markovian effects (they are all 
related to the phase of the evolving qubit) and can be 
safely neglected in this study. The interaction Hamilto- 
nians, instead, have a small but crucial difference, arising 
from the different trapping potentials of the impurities: 
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where no is the cond ensate density, i 7 is the quantization 
volume, Ek = yj ek(ek + 2rio.9s) is the Bogoliubov 
dispersion relation, ek = 7i 2 fc 2 /(2ms) is the dispersion 
relation of a non-interacting gas with k — |k| and Cfc 
is the Bogoliubov excitation operator. Operator n is 
the number operator of the impurities: For model I 
we assume that there is exactly one impurity atom in 
the double well system and therefore hn = |(1 + cr z ) 
and fiL = ^(1 — cr z ): where a z = |Z) (l\ — \r) (r|. The 
two wells are spatially separated by distance L so that 
x^ = X£ — L. For model II we also assume one impurity 
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in the lattice site. The atom has one internal state \g) 
which decouples from the environment and one |e) which 
does not and therefore h = \e) (e\. 

We note that H^g del 1 effectively describes two spatially 
separated qubits of Model II, albeit with a restricted 
state space {\eg) , \ge}}. Interestingly these states span 
the so called subdccohercnt state \eg) + \ge), which 
is very robust a gain st dcphasing noise induced by the 
environment (23l. |24| . Therefore we can expect the qubit 
architecture of model I to be less affected by noise than 
model II. 



III. REDUCED DYNAMICS AND 
INFORMATION FLOW 

The reduced dynamics of both models can be solved 
analytically 0, [25| . Each qubits dephases under the ef- 
fect of the ultracold gas, i.e., the diagonal elements of 
the qubit remain constant while the off-diagonals decay 
as p i(t) = e- r W+ ie (*Voi(0). The phase 9(t) has no 
effect on the information flow and therefore we do not 
consider it in this work. Instead we focus on the deco- 
herence function T(t). When T'(t) > information flows 
from the system to the environment and if there is an 
interval where T'(t) < then the flow of information 
is temporarily reversed. We associate this reversal with 
non-Markovian effects, adopting the proposal of Brcuer 
et al. as our definition for non-Markovianity Indeed, 
this measure of non-Markovianity, applied to a dcphasing 
model such as the two models considered here is 



,14 (a) ID-environment 
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Recall that this measure captures the maximal amount 
of information that can flow back from the environment 
to the system. The decoherence functions for the two 
physical systems considered here are 
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where a is the variance parameter. Interestingly the de- 
coherence factor of Model I has exactly the structure of 
the decoherence factor of two qubits of Model II with 
spatial separation L in a subdecoherent state p2j . as we 
anticipated in the previous Section. Therefore we can ex- 
pect some " coherence trapping" in model I that we would 
not observe in Model II. 



A. Dynamics of the decoherence factor 

We plot the decoherence factors of the two models in 
Fig. [5J using the same values of parameters as in Rcf. 
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FIG. 2. Decoherence functions r(t) Modcl 1 (solid line) and 
r ^ Model ii ( dashed line ) for ( a ) one-dimensional, (b) two- 
dimensional and (c) three-dimensional environment for as = 
0.25a Rb (black lines) and os = Oflb (gray lines). 



[Ta | and going to the limit of a continuum of modes, 
^ k (2tt) d J dk, where D is the dimension of the 
BEC. As anticipated, the qubit of Model I is much more 
robust against decoherence than the qubit of Model II. 
This difference is most striking in the case of a quasi- 
1D environment, where r(£) Model 1 saturates quickly to a 
small value, while r(£) Modcl 11 increases without bound 
over all considered time-scales. Hence Model I is al- 
most unaffected by the environmental noise while Model 
II loses all coherence, and all initial states converge to- 
wards the maximally mixed state. The difference in the 
dynamics of the two models is less drastic when the en- 
vironment is quasi-2D or 3D, where both decoherence 
factors saturate to a stationary value, although also in 
these two cases Model I is more robust against the noise 
than Model II. 

Furthermore, the decoherence factor of Model II is mono- 
tonic in the case of ID and 2D background gases, im- 
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plying that the flow of information is always from to the 
qubit to the environment and the dynamics is Markovian. 
This is at variance with the decoherence of Model I where 
we find both Markovian and non-Markovian dynamics in 
the case of ID and 2D background gases, depending on 
the value of the scattering length of the background gas 
[l5| . In the case of a 3D background gas, shown in Fig. 
[IJc), both decoherence factors are non-monotonic for a 
large enough value of the scattering length. This signals 
non-Markovian effects. In the next Section we quantify 
these using the measure of Eq. §5§ . 

B. Non-Markovianity 

In all cases we have considered we only see one pe- 
riod of information backflow. This permits the use of a 
slightly modified non-Markovianity measure, which cap- 
tures the maximal fraction of information that can flow 
back from the environment to the system after an ini- 
tial period of information flowing from the system to the 
environment [ToT ]. In Fig. [3] we show this modified non- 
Markovianity measure against the scattering length of 
the background gas in the case when the qubits are im- 
mersed in a 3D background gas. We observe a crossover 
from Markovian to non-Markovian dynamics for both 
models, although with different values of the crossover 
point: for a range of scattering lenghts Model I decoheres 
in a non-Markovian way, while the dynamics of Model II 
is Markovian. We also consider the effect of temperature 
on both systems. For a small temperature T = 10 nK we 
still observe a Markovian to non-Markovian crossover for 
Model I, although the crossover point has been shifted 
slightly to a larger scattering length. Physically this 
means that in order to induce non-Markovian dynam- 
ics, the boson-boson interaction has to be stronger to 
overcome the detrimental effects of thermal fluctuations. 
For Model II the measure is zero (the dynamics is Marko- 
vian) for all the considered values of the scattering length, 
demonstrating the non-Markovianity of Model II is very 
fragile against thermal effects in the environment. We ex- 
plain the differences in the dynamics of the two models 
in the following Section. 

IV. SPECTRAL DENSITY FUNCTIONS 

The spectral density function J(u>) — J3k I.9k| 2 i5(w — 
£k) characterising the dephasing dynamics of an open 
quantum system is determined by the coupling con- 
stants <7k of the effective interaction Hamiltonian Hab = 
°z Ylk 9kb^+H.c. The spectrum of each qubit model con- 
sidered in this article provides a framework for explain- 
ing the notable differences in their dynamics, namely the 
Markovian to non-Markovian crossover and the increased 
robustness against environmental noise of Model I com- 
pared to Model II. The spectral density, for small fre- 
quencies, shows a power-law behaviour J(u>) a w s . In 
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FIG. 3. Non-Markovianity measure for model I in a zero-T 
reservoir (solid black line) and a T = 10 nK reservoir (solid 
gray line) and model II in a zero-T reservoir (dashed black 
line). 
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FIG. 4. Spectral density functions J(lj) for Model I (solid 
lines in all figures) and Model II (dashed lines in all figures) 
in a (a) one-dimensional, (b) two-dimensional and (c) three- 
dimensional environment. Left hand side figures show the full 
spectrum, and the figures on the right show the low-frequency 
contribution. In the latter we show the spectrum for a weakly 
interacting background gas with as = 10 _3 a_Rb (black lines) 
and for a BEC with as = am (gray lines). 



the following we call parameter s > the Ohmicity pa- 
rameter since it determines whether the spectral density 
is Ohmic with s = 1, sub-Ohmic with s < 1 or super- 
Ohmic with s > 1. 

We conjectured in Ref. that the existence of the 

crossover from Markovian to non-Markovian dynamics 
is closely related to which class of the Ohmic spec- 
tral densities (sub-Ohmic, Ohmic or super-Ohmic) the 
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spectrum belongs to. In Ref. [21j we further quan- 
tified this claim, presenting a necessary condition for 
non-Markovian dephasing dynamics: Non-Markovian dy- 
namics can only appear if the spectrum of the envi- 
ronment is super-Ohmic. More specifically, we showed 
that for a very general dephasing model introduced by 
Palma, Suominen and Eckcrt, with a spectrum of the 
form J(uj) = ri u>^ s uj s exp{— w/oj c }, where 77 is a dimen- 
sionless coupling constant that we set to unity and ui p h 
is a phononic reference frequency introduced to keep the 
dimension of the spectrum equal for all values of s) 0, 
non-Markovian effects take place only if s > s cr it = 2, 
i.e., the spectrum is strongly super-Ohmic. For the phys- 
ical models considered here the form of the spectral den- 
sity function is more complicated, but we show that the 
main results hold also for these two systems. 



A. Spectra at low frequencies 

We plot the spectral densities J(w) for the two mod- 
els in Fig. HI focusing on the low-frequency part of the 
spectum in the right hand side column. We observed in 
Ref. I5[ that for the ultracold environments the effective 
Ohmicity parameter depends on the scattering length of 
the environmental bosons, s = s(ob) and on the dimen- 
sionality of the BEC environment. Changing these two 
allows transitions from one Ohmic class to another. Here 
we confirm that for both models and for all three di- 
mensions increasing the scattering length increases the 
effective value of s for small frequencies u. This effect 
is especially pronounced in the case of Model II. When 
the environment is a quasi-lD free gas (clb /am <C 1) the 
low-frequency part of the spectrum is sub-Ohmic but as 
the scattering length is increased to 0,3 / a-Rb ~ 1 ; the 
spectrum approaches an Ohmic form. However, even 
with further increase in the strength of the boson-boson 
coupling in the environment the spectrum does not be- 
come super-Ohmic and indeed we never observe non- 
Markovian effects in Model II in the ID regime. In the 
case of a quasi-2D environment the spectrum of the envi- 
ronment changes from sub-Ohmic to super-Ohmic when 
the scattering length value is increased, but even in the 
more strongly interacting case the spectrum is not suffi- 
ciently super-Ohmic to trigger non-Markovian effects. In- 
stead in the 3D case the free gas has an Ohmic spectrum 
which turns super-Ohmic as the interaction between the 
bosons is turned on and increased. In this case the spec- 
trum can become so super-Ohmic that we observe non- 
Markovian effects. 

Model I is already naturally more Ohmic than Model II 
in the sense that when the ultracold gas is essentially 
non-interacting both the quasi-2D and the 3D environ- 
ments have a super-Ohmic spectrum (the quasi ID en- 
vironment is roughly Ohmic) and the spectra become 
super-Ohmic in all dimensions when the scattering length 
is increased. Critically, in all dimensions the spectrum 
becomes super-Ohmic enough to create non-Markovian 



effects in the dynamics of the double well qubit for some 
critical value of the scattering length, and therefore we 
observe a crossover from Markovian to non-Markovian 
dynamics for Model I in all three dimensions. 



B. Full spectrum 

Changes in the scattering length have a crucial effect 
on the spectral density for small frequencies but they 
are almost negligible when looking at the full spectra. 
Instead the full spectra, shown in the left-hand-side col- 
umn of Fig. 21 exhibits another very interesting difference 
between the two models, namely that the spectrum of 
Model I oscillates as a function of tu. Moreover, the spec- 
tral density function vanishes for some specific values of 
lu in the quasi-lD case, implying that some modes of the 
environment are completely decoupled from the qubit. 
We find that the larger is the separation between the two 
wells the more roots the spectral density has, i.e., more 
modes decouple from the qubit. We observed numeri- 
cally that increasing the well separation also increases 
the " coherence trapping" , leading to higher stationary 
values of the off-diagonal elements of the qubit density 
matrix. The higher is the dimension of the environment, 
the smaller are the deviations of the spectrum of Model 
I compared to that of Model II. This phenomenon is also 
reflected in the differences in the decoherence factors of 
the two models; the differences arc most pronounced in 
the case of a quasi-ID environment and in the higher 
dimensions the decoherence factors are more similar in 
both value and dynamical behaviour. 



V. DISCUSSION AND CONCLUSIONS 

We have studied the non-Markovian dynamics of two 
physically different realisations of a qubit interacting 
with a BEC environment. We discovered that the 
qubit architecture of Model I is much more sensitive to 
non-Markovian effects than the one used in Model II. 
This statement applies in three ways: (i) Model I has 
non-Markovian dynamics in all three dimensions of the 
BEC, unlike Model II, which is Markovian in ID and 
2D environemnts: (ii) it has larger values of the non- 
Markovianity measure in the cases when both qubits have 
non-Markovian dynamics, i.e., it regains previously lost 
information more easily: (iii) it is more robust against 
thermal fluctuations. We also discovered that the two 
qubit architectures can have extremely different sensitiv- 
ity to environmental noise, especially in the case when 
the environment is effectively one-dimensional. This re- 
sult highlights the importance of choosing suitable qubit 
systems when designing quantum simulations and quan- 
tum probe systems. 

We also explored the connection between the form of the 
spectral density function and the ensuing qubit dynam- 
ics. While the form of the full spectral density function 
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dictates the general dynamics of the qubit, only a very 
small part of it controls the Markovian to non-Markovian 
crossover. We noted the importance of the low-frequency 
part of the spectrum in the emergence of non-Markovian 
phenomena in Ref. [2l| in the case of a simple dephas- 
ing model, and the study presented in this article con- 
firms that the statement holds also for more complex 
systems. It is nonetheless striking to notice the over- 
whelming importance of the low-frequency modes. The 
spectral density function of model I in a quasi-lD BEC 
has a very rich structure over the frequency range of the 
order of the cut-off frequency er -1 , yet this has no effect 
on the crossover of the qubit dynamics from Markovian 
to non-Markovian. The crossover is fully controlled by 
the behaviour of the spectral density function for low fre- 
quencies oj <C c" 1 , specifically whether the spectrum is 
quadratically increasing or not. It is worth stressing that 
this connection seems to be quite specific to pure dephas- 
ing noise. In Ref. |26[ the authors studied a dissipative 
model and found a direct connection between roots of 
the spectral density function and non-Markovian dissi- 
pationlcss dynamics; here we discovered that in the pure 
dephasing qubit model the roots of the spectra do not at 
all affect the (non-)Markovianity of the dynamics. 



In conclusion, our results have twofold importance. On 
one hand they illustrate in a very clear way how the con- 
nection between non-Markovianity (or memory effects) 
and structured environments generally has to be taken 
with great care. On the other hand, and most impor- 
tantly, they warn us of the misuse of the term "non- 
Markovian environment" . In our study the environment 
is exactly the same for the two models, and in both cases 
it is interacting with a qubit probe. However, under cer- 
tain conditions, perfectly identical environments induce 
Markovian dynamics on one qubit and non-Markovian 
dynamics on another. 



ACKNOWLEDGMENTS 

We acknowledge financial support from EPSRC 
(EP/J016349/1), the Finnish Cultural Foundation (Sci- 
ence Workshop on Entanglement) and the Emil Aaltonen 
foundation (Non-Markovian Quantum Information) and 
Magnus Ehrnroth foundation. We thank Gabriele De 
Chiara, Dieter Jaksch, Stephen Clark and Tomi Johnson 
for helpful discussions. 



[1] H. P. Breuer and F. Petruccione, The Theory of Open 
Quantum Systems, (Oxford Univ. Press, 2007). 

[2] U. Weiss, Quantum Dissipative Systems (World Scien- 
tific, Singapore, 1999). 

[3] G. Lindblad, Commun. Math. Phys. 48, 119 (1976). 

[4] M. M. Wolf and J. I. Cirac, Commun. Math. Phys. 279, 
147 (2008). 

[5] H.-P. Breuer, E.-M. Laine and J. Piilo, Phys. Rev. Lett. 

103, 210401 (2009). 
[6] A. Rivas, S. F. Huelga and M. B. Plenio, 

Phys. Rev. Lett. 105, 050403 (2010). 
[7] X. M. Lu, X. Wang and C. P. Sun, Phys. Rev. A. 82, 

042103 (2010). 

[8] R. Vasile, S. Maniscalco, M. G. A. Paris, H.-P. Breuer, 
and J. Piilo, Phys. Rev. A 84, 052118 (2011). 

[9] S. Luo, S. Fu, and H. Song, Phys. Rev. A 86 044101 
(2012). 

[10] P. Haikka, J. D. Cresser, and S. Maniscalco, Phys. Rev. 

A 83, 012112 (2011) 
[11] L. Mazzola, E.-M. Laine, H.-P. Breuer, S. Maniscalco, 

and J. Piilo, Phys. Rev. A 81, 062120 (2010); E.-M. 

Laine, K. Luoma, and J. Piilo, Jour. Phys B 45, 154004 

(2012). 

[12] E.-M. Laine, J. Piilo, and H.-P. Breuer, Europhys. Lett. 

92, 60010 (2010); A. Smirne, D. Brivio, S. Cialdi, B. 

Vacchini, and M. G. A. Paris, Phys. Rev. A 84, 032112 

(2011); C.-F. Li, J.-S. Tang, Y.-L. Li, and G.-C. Guo, 

Phys. Rev. A 83, 064102 (2011). 
[13] P. Haikka, J. Goold, S. McEndoo, F. Plastina and S. 

Maniscalco, Phys. Rev. A 85, 060101 (2012). 
[14] B.-H. Liu et at, Nature Physics 7, 931 (2011). 
[15] P. Haikka, S. McEndoo, G. De Chiara, G. M. Palma and 

S. Maniscalco, Phys. Rev. A 84 031602(R) (2011). 



[16] I. Bloch, J. Dalibard, and W. Zwerger, Rev. Mod. Phys. 
80, 885 (2008); C. Chin, R. Grimm, P. Julienne, and E. 
Tiesinga, Rev. Mod. Phys. 82, 1225 (2010). 

[17] See, for example, T. Tufarelli, A. Ferraro, M. S. Kim, 
and S. Bose, Phys. Rev. A 85, 032334 (2012); J. Goold, 
T. Fogarty, N. Lo Gullo, M. Paternostro, and Th. Busch, 
Phys. Rev. A 84, 063632 (2011); S. Campbell, M. Pa- 
ternostro, S. Bose, and M. S. Kim, Phys. Rev. A 81, 
050301 (2010), S. Fishman, G. De Chiara, T. Calarco, 
and G. Morigi, Phys. Rev. B 77, 064111(2008) and A. 
Recati, P. O. Fedichev, W. Zwerger, J. von Delft, and P. 
Zoller, Phys. Rev. Lett. 94, 040404 (2005). 

[18] A. Klein, M. Bruderer, S. Clark and D. Jaksch, New J. 
Phys. 9, 411(2007). 

[19] M. Bruderer, A. Klein, S. Clark and D. Jaksch, Phys. 
Rev. A 76 011605(R) (2007). 

[20] J.-G. Li, J. Zou and B. Shao, Phys. Rev. A 81, 062124 
(2010). 

[21] P. Haikka, T. H. Johnson and S. Maniscalco, 

larXiv: 1203.64691 
[22] M. A. Cirone, G. De Chiara, G. M. Palma and A. Recati, 

N. Jour. Phys. 11,103055 (2009). 
[23] G. M. Palma, K.-A. Suominen and A. K. Ekert, Proc. 

Roy. Soc. Lond. A 452 567 (1996). 
[24] R. Doll, M. Wubs, P. Hanggi and S. Kohler, Europhys. 

Lett. 76, 547 (2006). 
[25] R. Doll, D. Zueco, M. Wubs, S. Kohler, and P. Hanggi, 

Chem. Phys. 347, 243 (2008). 
[26] W.-M. Zhang, P.-Y. Lo, H.-N. Xiong, M. W.-Y. Tu and 

F. Nori, Phys. Rev. Lett. 109, 170402 (2012). 



